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ABSTRACT 


Synge's  formulation  for  Hamilton's  theory  of  geometric  optics 
is  used  to  specify  the  rays  for  a  variety  of  problems.  Based  upon 
the  analogy  between  ray  and  particle  trajectories,  an  attempt  is 
made  to  treat  the  problem  of  beyond- the-horizon  atmospheric  propa¬ 
gation.  It  is  shown  that  the  observed  propagation  law  (e~ks)  can 
be  motivated  from  the  concept  that  the  received  signal  is  the  sum 
of  the  Hamilton  rays  linking  the  transmitter  and  receiver.  This 
implies  three  differences  with  respect  to  the  Booker-Gordon  theory 
for  the  everpresent  scatter  signal:  a  change  in  the  correlation 
function  from  exponential  to  Gaussian;  adoption  of  millimeters  as 
the  size  of  the  characteristic  fluctuations;  and  an  extension  of 
the  critical  volume  to  include  manyfold  scatter  in  the  tropospheric 
path  from  the  optical  horizon  to  the  receiver. 
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SECTION  I 


INTRODUCTION 

The  ray  treatment  also  called  geometric  optics,  is  considered 
to  be  appropriate  when  the  wavelength  is  small  compared  to  the 
characteristic  length  in  slowly  varying  media  comprising  the  optical 
instrument,  yet  large  compared  to  molecular  dimensions.  Within  the 
realm  of  geometric  optics,  there  are  several  possible  starting  points. 

For  example,  a  generalized  form  of  the  Sommerfeld-Runge  relation 
(equivalent  to  Snell's  formulation)  or  forms  derivable  from  the 

Hamilton-Fermat  principle  (the  Hamilton- Jacobi  differential-equation 
formulation  or  Hamilton's  canonical  equation)  can  be  employed. 

Until  recently,  the  first  system  suffered  because  it  was  very 
difficult  to  formulate  Snell's  law  for  the  most  general  case  of 
inhomogeneous  and  anisotropic  media.  However,  this  situation  has 
apparently  been  alleviated  since  Poeverlein  [Ref.  l]  has  recently 
shown  that  the  Sommerfeld-Runge  "Law"  (V  x  K  -  0)  together  with 
F(r,  k)  =  0  provide  for  the  rays  under  a  wide  range  of  conditions, 
including  anisotropy  and  inhomogeneity.  Absorption  and  lateral- 
intensity  variations  can  be  included  by  introducing  complex  pro¬ 
pagation  vectors.  Further  generalization  to  four  dimensions  en¬ 
ables  treatment  of  modulated  waves  and  time-varying  media.  Once 
formulated,  calculation  of  the  rays  should  be  straight  forward. 

•k 

Numbers  in  bracket  designate  References  listed  at  end  of  report. 
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However,  conservation  of  the  tangential  components  across  discon¬ 
tinuities  appears  to  be  afflicted  with  the  same  limitation  to  strati¬ 
fied  media  that  plagued  the  simpler  theory. 

On  the  other  hand,  the  Hamiltonian  system  permits  almost  routine 
formulation  of  the  most  general  problem,  and  there  is  no  limitation 
to  stratification  or  to  constant  B-field.  For  these  reasons,  and 
since  Hamilton's  original  work  was  performed  approximately  a  cen¬ 
tury  ago,  it  is  difficult  to  account  completely  for  the  fact  that 
Hamilton’s  optics  is  not  universally  taught.  Perhaps  one  reason 
is  that  there  is  no  motivation  in  optics  comparable  to  the  role 
of  Hamilton's  discipline  as  a  prerequisite  for  quantum  mechanics. 
However,  the  power  and  elegance  of  Hamilton's  optics  alone  are 
sufficient  justification  for  more  widespread  use. 

Since  the  system  usually  leads  to  a  set  of  nonlinear  second- 
order  differential  equations,  a  computer  is  the  only  convenient  means 
for  computation  of  the  rays;  this,  however,  is  true  of  any  discipline, 
provided  the  medium  is  sufficiently  arbitrary.  To  date,  the  large 
digital  computer  can  best  handle  the  anisotropic  ray  computation. 

Wong  [Ref.  2]  has  employed  the  analog  computer  for  the  troposphere 

and  is  currently  synthesizing  a  program  for  using  this  device  for 
the  more  difficult  (inhomogeneous  anisotropic)  ionosphere.  Un¬ 
fortunately,  it  appears  that  without  extensive  modification  the 
analog  computer  does  not  carry  enough  significant  figures  to  give 
reproducible  results. 
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For  purposes  of  this  report,  it  is  believed  that  that  Hamil¬ 
tonian  system  will  provide  the  best  means  of  attacking  new  problems 
since  the  analogous  disciplines  in  mechanics,  statistical  mechanics, 
and  quantum  mechanics  can  be  used  as  a  guide. 
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SECTION  II 


RAY  TREATMENT 

HAMILTONIAN  OPTICS  FOR  CARTESIAN  COORDINATE  SYSTEMS 

For  Cartesian  coordinate  systems  the  Hamiltonian  system  as  given 

by  Synge  [Ref.  3]  is  particularly  simple  and  direct.  Consider  a  fixed 

(no  time- dependence)  optical  instrument  or  an  aggregate  of  media  whose 

phase  refractive  index  is  a  function  of  position  for  the  inhomogeneous 

case  and  a  function  of  direction  for  the  anisotropic  case.  Thus 

v  =  v  >  where  is  the  direction  cosine.  In  addition, 

a  a  =1,  where  the  repeated  index  implies  the  summation  convention. 

K.  K. 


The  Euler-Lagrangian  condition  that  the  ray  path  with  fixed 
end-points  must  be  an  extremum,  as  obtained  from  the  variational 
principle  6  J  (l/v)ds  =  0  is  expressed  by 

±  jL  (1)  _  _fL  (I)  o  . 

ds  d ctk  \  v  /  dxk  \  v  / 


Because  of  the  auxiliary  condition  =  it  is  possible 


to  decide  further  that  the  Langrangian  (1/v)  shall  be  homogeneous 

a 


of  degree- plus -one  in  the  Then  obviously, 


dor. 


(1/v)  is 


of  degree-zero  in  the 


This  expression  defines  the  canonical 


variable  of  slowness, 
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Figure  1  gives  the  relation  between  the  ray  surface  or  indica- 
trix  r  =  v  ^x^>  anc*  the  basic  wave  surface  or  figuratrix 
^  =  0.  These  are  the  reciprocal  surfaces  of  Hamilton 
and  are  related  by  PS  •  PR  =  PS  .  PQ  cos  6  =  1 


In  principle,  the  ray  paths  can  be  derived  from  either  of  the  two 
surfaces;  however,  one  is  usually  given  the  phase-refractive  index 
of  the  medium.  Since  the  figuratrix  or  ray  Hamiltonian  is  given 
in  terms  of  the  phase  refractive  index,  it  is  more  convenient  to 
employ  the  ray  Hamiltonian  to  derive  the  rays.  Indeed,  vanishing 
of  the  ray  Hamiltonian  and  the  system  of  canonically  conjugate 
equations , 

=  =  _  ifk  (2) 

da jj  du  '  du 
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provides  a  sufficiency  condition  for  determining  the  rays.  The 


I 


parameter  u  is  related  to  the  time  t  through  the  degree 


of  homogeneity  of  (ft  +  1)  in  ,  or  more  explicitly, 


du  . 


in  the  case  of  discontinuous  media,  by  the  jump  conditions 


=  0  • 


(3) 


In  addition,  there  is  the  alternative  of  replacing  the  dis¬ 
continuity  with  a  thin  continuous  layer  that  provides  the  ap¬ 
propriate  ray  trajectory  consistent  with  Equations  (2)  and  (3). 

GENERATING  FUNCTIONS  OF  THE  RAY  HAMILTONIAN 

The  generating  function  (Synge)  appropriate  to  both  inhomo 
geneous  and  anisotropic  media,  but  restricted  to  Cartesian  co¬ 
ordinate  systems  is  given  by 


n  (xk  -  °k>  = 


V  gkgk 


1  =  0, 


(4) 


f  (xk,  al/a2,  al/°l) 
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where  c 


1  throughout.  In  this  case  Q  +  1  is  homogeneous 


of  degree-plus-one  in  the  ,  so  that,  further,  u  =  t.  To  be 
able  to  employ  orthogonal  systems  other  than  the  Cartesian  system, 
it  is  acceptable  to  set  up  the  ray  Hamiltonian  in  the  above  form 
and  then  to  transform  coordinates  in  the  final  system  of  differen¬ 
tial  equations. 

In  the  special  case  of  isotropic  media,  the  ray  Hamiltonian 
for  all  orthogonal  coordinate  systems  is  given  by  (Synge)* 


glk  <Tk  -  n2  (Xk)  =  0  , 


C  5) 


where  n  is  the  phase-refractive  index  and  g  is  the  con- 


travariant  metric  fundamental  tensor  defined  by  the  familiar 
relations 


ik 


minor  of 


g 


^Private  correspondence  with  the  author. 


The  Apple ton-Hartree-Goubauian  formulation  of  the  refractive 


index  for  the  ionosphere  is  given  by 


If  the  ion  contribution,  the  collision  frequency,  and  the  Lorentzian 
polarization  term  can  be  ignored,  the  refractive  index  reduces  to 


n 


2 


1 


1 


,  2 

r  2 

yT 

yT 

'  2(1  + a)  1 

4(1  + a'" 

(7) 


and  the  ray  Hamiltonian  describes  an  e1  psoid  of  revolution  whose 

axis  of  symmetry  is  along  the  terres  .al  magnetic-field  lines.  If 

the  further  restriction  is  made  tl  „  radio  frequencies  are  in  excess 

of  median  MUF^  (maximum  usable  frequency  for  path  s  ),  then 
2  2 

^  ^  and  the  expression  for  the  index  reduces  to 


n2  .  l 


4  7T  N_ 


m  co 


=  1  - 


cor 


(8) 
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In  this  case,  the  ellipsoidal  surface  defining  the  ray  Hamiltonian 
degenerates  to  a  sphere  appropriate  to  isotropic  media. 

NOTE 


It  is  now  customary  to  ignore  the  Lorentzian  polarization  term 

(4/3  ttP)  in  the  contribution  of  free  electrons  to  the  refractive 

index  of  the  ionosphere.  It  is  included  in  the  general  expression, 

Equation  (7),  since  it  is  evident  that  the  term  is  pertinent  to  the 

ion  contribution.  The  Lorentzian  term  takes  into  account  the  self- 

field  of  the  specie  as  polarized  by  the  external  field.  The  electron, 

-13 

a  monopole  with  a  dimension  of  the  order  of  10  cm,  moves  with  ease 
to  minimize  the  distortion.  In  the  case  of  the  ion  (a  complex  ion 
such  as  0+,  and  not  H+)the  charge  separation  (of  the  order  of 
10"  cm  )  and  the  relative  immobility  of  the  ion  permit  the  specie 
to  be  polarized  to  such  a  degree  that  the  self-field  is  an  important 
contribution. 

THE  TENSOR  FORM  OF  THE  VARIATIONAL  PRINCIPLE 

The  presentation  of  Hamilton's  system  as  formulated  by  Synge 
began  with  an  expression  of  the  variational  principle  which,  while 
obviously  applicable  to  inhomogeneous  and  anisotropic  media,  was 
limited  to  Cartesian  coordinate  systems.  For  the  most  general 
media  it  was  possible  to  take  advantage  of  preferred  coordinate 
systems  or  ignorable  coordinates  only  by  finally  transforming  the 
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system  of  differential  equations.  Such  a  scheme  does  not  take  full 
advantage  of  the  ignorable  coordinate.  Haselgrove  ^Ref .  4]  has  demon¬ 
strated  that  a  start  can  be  made  with  a  tensor  expression  of  the  var¬ 
iational  principle  which  is  appropriate  to  all  orthogonal  coordinate 
systems  and,  of  course,  to  both  anisotropic  and  inhomogeneous  media. 

By  replacing  total  and  partial  derivatives  by  terms  involving  the 
covariant  derivative.  Equation  (1)  becomes 


where  the  reciprocal  surfaces  of  Hamilton  are  F  =  1  and  G  =1; 

j  H i,m  |  denotes  the  familiar  Christoffel  symbol  of  the  second  kind; 
and  x^  and  u^  donote  the  position  and  velocity,  respectively. 

The  double  notation^  r  =  v,  Q  =  0;  F  =  1,  G  =  1^  has  been 
retained  because  the  surface  can  be  expressed  either  in  terms  of  the 
slowness  vector  or  its  inverse,  the  velocity  vector.  Hamilton’s 
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canonical  forms  can  then  be  written  as 


(90 

<9u* 


^  =  gij  1SL  _  „ij  „km  u£  (9G  /  (jg£i  _ 
dt  dxi  dum  \dxk  ) 


(10) 


Equations  (10),  together  with  the  basic  wave  surface  G  =  1  define  the 
rays  for  any  coordinate  system.  For  the  spherically  symmetric  case 
the  six  Equations  (1)  reduce  (Haselgrove)  to 


dvr 

dt 


( A  is.  +  vo  *o 

V  n  dt  dt 


+  sin  0  v 


d<l>\ 
dt  ) 


1  (l_  dn_ 
r  \  n  dd 


r  cos  6  v 


<t> 


d$\ 
dt  ) 


dv^  =  1  /  1_  dn_ 

dt  r  sin  6  \  n  <9<I> 


sin  d 


dr 

dt 


r  cos  0  v 


<J>  d 0 
dt 


(11) 
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In  spite  of  their  tensor  form,  Equations  (9)  and  (10)  do  not  define 
a  Lorentz  covariant  theory.  This  is  clear  because  the  frequency 
enters  only  in  determination  of  the  refractive  index  n  and  not 
as  a  variable  having  the  status  of  x_.  ,  u^  .  A  four-dimensional 
discipline,  allegedly  appropriate  to  a  wider  variety  of  problems, 
has  been  formulated  by  Synge  (see  pages  25  to  28).  The  Haselgrove 
discipline  could  be  employed  for  some  problems  by  taking  the  re¬ 
fractive  index  at  the  appropriate  value  of  the  Doppler  frequency  as 
given  by  relatively. 


THE  SPHERICALLY- SYMMETRIC  INHOMOGENEOUS  ISOTROPIC  CASE 

In  Hamiltonian  optics,  as  in  Hamiltonian  mechanics,  the  co¬ 
ordinate  system  is  chosen  to  take  advantage  of  the  existence  of 
ignorable  coordinates,  since  the  latter  lead  immediately  to  con¬ 
stants  of  the  ray  or  particle  trajectory.  Invariably,  the  canonical 
conjugate  of  an  ignorable  coordinate  is  a  constant  of  the  motion. 

In.  the  isotropic  case  the  ray  Hamiltonian  is  given  by  Equation  (5), 
and  the  coordinate  system  is  chosen  so  that  g  and  n  involve 

1.  K. 

the  smallest  number  of  coordinates.  The  selected  differential  line 
element  is  then  ds2  =  dS*  d§^  =  dr2  +  r2d02  +  r2  sin  0  d§2  , 

and  the  ray  Hamiltonian  takes  the  form  (again  c  =  1) 


12 


n  =  -  H2  (Xk>  =  0 


n  = 


1 

r2  sin2  6 


n2  (t)  =  0  . 


(12) 


The  origin  is,  of  course,  coincident  with  the  center  of  symmetry. 

It  is  worth  repeating  that  vanishing  of  the  ray  Hamiltonian  and  the 
canonical  equations  is  sufficient  to  uniquely  define  the  ray  tra¬ 
jectories.  The  canonical  Equations  (2)  give 


dr 

du 


=  2a 


r 


d  d  ^6_ 
du  "  r2 


d(J)  2(7<f)  2C1 

r2  sin^  6  r2  sin^  0 


d  a 


r 


du 


2a0  2 

r3 


+ 


2<r<t)  2 
sin2  6 


+ 


2n 


dn 

dr 


2ff(j)2  COS  e 
r2  sin^0 


du  =  0;  oq  =  C2  . 


(13) 
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canonically  conjugated  to  the 


As  anticipated,  the  coordinate 
ignorable  coordinate  $  is  a  constant  of  the  ray  trajectory.  In 
addition,  inspection  of  the  homogeneity  of 

11  (r,  6,  oT,  gq,  (7$)  =  0  in  o ^ 

indicates  that  the  phase  parameter  u  =  t/2.  The  number  of  equa¬ 

tions  can  be  reduced  by  eliminating  the  o^,  although  this  form  is 

not  necessarily  more  advantageous  for  computer  use. 
c  ...  .p,  dr  „  d0  d$  , 

Settins  p  =  Q  =  to'  R  =  dT!  Produces 


dP 

du 


rQ2  + 


4  Ci 


r2  sin2  0 


4  n 


dn 

dr 


d_ 

du 


(r2Q) 


4C^2  cos  0 
r2  sin^  9 


2C1 

r2  sin2  9 


(14) 


For  application  to  a  spherical  electron  cloud  with  radial  symmetry, 
solution  of  Equation  (14)  yields  a  three-dimensional  picture  of  the 
rays.  The  simultaneous  solution  of  properly  selected  rays,  together 
with  the  equation  of  the  terrestrial  surface,  can  also  provide  an 
estimate  of  the  ray  pattern  or  relative  intensity  of  radiation  on 
the  ground.  The  effect  of  absorption  along  the  ray  trajectory  in 
the  ionized  cloud  can  also  be  included. 
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If  infinite  skew-rays  are  ignored  the  problem  of  defining  ray 
paths  in  a  spherical ly- symmetric  medium  (ion  cloud)  can  also  be  solved 


in  two  dimensions  by  taking  a  planar  cut  through  the  cloud.  The 
plane  is  defined  by  two  coincident  lines:  the  initiating  ray,  and  a 
line  between  the  point  of  ray  initiation  and  the  center  of  symmetry. 
Polar  coordinates  are  again  selected  to  take  advantage  of  the  exis¬ 
tence  of  an  ignorable  coordinate,  but  it  is  to  be  emphasized  that  only 
the  0  coordinate  has  meaning  in  a  particular  plane.  The  ray  Hamil¬ 
tonian  becomes 


or2  +  ~Y  °o 2  -  "2  «  =  0 ' 


(15) 


r 


and  the  Hamiltonian  canonical  equations  give  the  four  equations 


d  0  2C3 

du  r2 


(16) 
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The  homogeneity  of  the  ray  Hamiltonian  in  has  not  changed,  so 

t  =  2u.  Substituting  for  the  and  \  immediately  produces 


But,  since 


£  (t)  *^32  £  (£)  -2*  (»2)-0. 

0  is  ignorable,  d/du  =  (dr/du)/  d/dr.  Integrating 

+  4C2r—  2  —  4n2  =  0,  (18) 


and  finally  for  the  ratio  of  d0/du  to  dr/du 


dO 

dr 


.-1 


(Cn2  r2 


1)“*  • 


(19) 


It  is  well  known  that  this  expression  can  also  be  obtained  from  the 
cylindrical  form  of  Snell's  Law  (nr  sin  0  =  k),  but  not  with  the 
facility  and  generality  characteristic  of  the  Hamilton  system. 

RAY  PATHS  IN  ISOTROPIC  MEDIA  WITH  SPHEROIDAL  SYMMETRY 

This  particular  example  is  concerned  with  inhomogeneous  media 
having  equiref ractive  index  surfaces  that  are  prolate  or  oblate 
spheroids  (ellipsoids  of  revolution).  At  first  glance  it  might 
appear  that  ellipsoidal  coordinates  are  the  most  advantageous;  how¬ 
ever,  for  an  ignorable  coordinate  to  exist  in  the  ray  Hamiltonian, 
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neither  the  g#1  nor  the  n  can  contain  the  coordinate.  To 
ik 

satisfy  this  criterion  the  spherical  coordinate  system  is  again 
selected,  and  the  ray  Hamiltonian  is  given  by 


A  =  glk  ai  ^  -  nz  (xk)  =  0 


2  1  2 
°r  +  2  °0  + 


-X  1  9  a^2  -  n2  (r,  0)  =  0 

r2  sin2  0 


(  20) 


The  functional  dependence  of  the  refractive  index  upon  r,  0  ($  being 
ignorable  by  rotational  symmetry)  is  (e.g.,  artificial  cloud) 


=  1  —  A  exp 


{- 


n2  =  1  —  A  exp  .j—  i 


(1  - 


(1  - 


cos  0) j ~ |  £or  t^e  opiate  spheroid,  and 
i  sin  0)  *|  -  \  for  the  prolate  spheroid.  Again, 

P  J  * 

vanishing  of  the  ray  Hamiltonian  and  Hamilton's  canonical  equations 
is  sufficient  to  determine  the  ray  paths  as  expressed  in  Equation  (21) 


du 


■ae  n  a1>  <9n 

— ^ —  +  2  z  z  +  2  n  _ 
sin2  6  ^ r 


dae 


du 


2  cos  Q 
r-2  sin^  0 


2n 


dn 

dO 


d0(j) 

— -  =  0  ;  oq  = 


C 
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du  =  2ar 


de  _  2<7(9 

du  ~~  r2 


d(I> 

du 


2g6> 

r2  sin2  6 


(21) 


Inserting  the  values  of  and  ,  and,  further,  letting 

P  =  dr/du,  Q  =  d0/du,  and  R  =  d§/du,  the  system  of  equations 
reduces  to 


dP 

du 


=  r 


4C" 


r3  sin2  6 


4n 


<9n 

dr 


4C2  cos  Q 

r2  sin30 


4  n 


(9n 

dd 


d$ 

du 


2  C 

r2  sin2  d 


(22) 


INHOMOGENEOUS  AND  ANISOTROPIC  RAYS  IN  MEDIA  WITH  NO  SYMMETRY 

The  system  of  differential  equations  that  can  yield  ray  paths 
for  any  continuous  medium  to  which  geometric  optics  is  appropriate 
and  as  given  by  Haselgrove,  was  specified  in  Equations  (9)  and  (10) 
for  any  coordinate  system  and  in  Equation  (11)  for  polar  coordinates. 
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Generally  speaking,  an  anisotropy  in  an  ionized  gas  is  a  result  of  the 
impressed  magnetic  field.  Moreover,  the  field  lines  are  parallel  over 
the  extent  of  the  plasma  of  interest.  It  is  quite  natural,  then,  for 
the  specification  of  the  ray  Hamiltonian  in  the  form  of  an  ellipsoid 
of  revolution  to  be  simplest  in  Cartesian  coordinates  where  one  of  the 
coordinates  (say  z)  is  parallel  to  the  field  vector.  The  ray  Hamil¬ 
tonian  is  then  given  simply  as 


ft  (xk  ,  (7k) 


b2  (*k> 


c2  (xk) 


o  . 


(23) 


Again,  vanishing  of  the  ray  Hamiltonian  and  Hamilton’s  canonical 
equations  provides  a  sufficiency  condition  for  determining  the  ray 
trajectories. 

Substituting  the  values  of  and  the  cannonical  equations 

lead  to  the  system  of  nonlinear  second-order  differential  equations 


^  =  fx  (xk)  p2  +  gx  (xk)  Q2  +  hi  (*k)  r2 


=  f2  (xk)  P2  +  g2  (xk)  Q2  +  h2  (xk)  r2 


=  f3  (xk)  P2  +  g3  (xk)  Q2  +  h3  (xk)  R2 


(24) 
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The  system  can  also  include  time-dependent  media.  In  this  case, 
there  is  an  added  term  in  Equation  (24)  of  the  form  P  d/du  (f^),  etc. 
An  unsuccessful  attempt  was  made  to  linearize  the  system  of  equations 
so  that  the  modified  Peano-Baker  matrizant  methods  could  be  applied, 
although  the  system  in  its  existing  form  is  quite  amenable  to  computer 
techniques. 

ABSORPTION,  PHASE,  AND  TIME 

In  the  ray  treatment  to  this  point,  absorption  of  energy  has 
not  been  considered.  To  include  absorption,  the  time  increment  of 
group  propagation  must  be  known,  i.e.,  time- reckoning  along  the  ray 
path.  This  introduces  an  added  element  of  complexity  because  the 
presence  of  the  phase  refractive  index  in  the  previous  equations 
implies  that  the  time- reckoning  involved  is  the  phase  (periodicity 
of  the  wave).  The  true  time  is  then  given  in  terms  of  the  phase 
propagation  time,  t,  as 


dn 


n  +  oj  - 


(25) 


n 


phase 


Haselgrove  conveniently  includes  the  phenomenon  of  absorption  by 
introducing  the  above  equation,  together  with 


dD 

dt 


n 


n 


D  , 


(26) 
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where  h  =  k  (x^,  u  )  is  the  absorption  coefficient,  t* (c)  is 
the  equivalent  path  length  in  vacuo,  and  D  is  the  absorption 
D  =  exp  {  -  J  H  cos  0ds  }.  The  phase  t  is  consistently  used 
throughout;  and  the  simultaneous  solution  of  Equations  (9),  (10), 
(25),  (26)  gives  all  the  required  parameters  of  propagation  in  the 
most  general  case  of  time- stationary  media. 

Parenthetically,  the  confusion  between  phase  t  and  time 
t1  extends  through  much  of  physics.  For  example,  the  t  used  in 
most  of  quantum  theory  is  phase  not  time.  As  motivation  for  this 
assertion,  first  quantization  confers  a  wave  property  to  the  particle 
(defined  here  as  any  simply  connected  intracule  or  any  extracule, 
always  of  nonvanishing  rest  mass).  This  is  illustrated  by  experi¬ 
ments  on  a  queue  of  electrons  impacting  on  a  microscopic  double  slit. 
Intuitively,  the  electron  passes  through  one  slit  or  the  other. 
However,  the  accompanying  phase  wave,  propagating  at  a  velocity  in 
excess  of  the  characteristic  velocity  c  permits  the  electron  to 
be  aware  of  the  existence  of  both  apertures  and,  via  this  phase 
wave,  to  interfere  with  itself. 

It  might  be  argued  that  quantum  mechanical  space  is  everywhere 
homogeneous  and  isotropic  [see  Equations  (28),  (29)].  The  temporal 
reckonings  in  phase  t  and  time  t1,  therefore,  differ  only  by  a 
constant  factor  that  cannot  affect  the  essential  phenomenology. 
However,  Eddington  [Ref.  5]  has  asserted  that  the  difference  is 
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indeed  significant.  Quantum  mechanics  is  a  particle  theory,  i.e., 
one  in  which  the  perturbation  between  the  object  system  and  the  en¬ 
vironment  is  linked  to  the  object  system,  the  environment  remaining 
essentially  invariant  (the  Hartree-Fock  self-consisten  field). 
Consequently,  it  is  conceptually  possible  that,  in  certain  cases, 
the  perturbation  of  this  so-called  unpolarizable  environment  might 
be  sufficiently  strong  to  make  the  postulate  of  isotropic  homo¬ 
geneous  space  rather  unrealistic  and,  therefore,  a  relatively  incon¬ 
venient  starting  point. 

In  other  words,  a  reduced  electron  harmonically  oscillating 
in  a  radial  potential  well  may  present  an  environment  quite  differ¬ 
ent  from  that  of  a  reduced  nucleon  also  harmonically  oscillating 
in  the  nucleus.  Consequently,  interpretation  of  the  quantum- mechan¬ 
ical  t  may  actually  depend  upon  the  physical  problem,  and  may  be¬ 
come  a  vital  consideration  in  nucleon  phenomenology. 

If  further  motivation  for  the  existing  state  of  affairs  is 
required,  historically,  the  most  vital  problem  solved  by  the  quantum 
mechanist  is  =  e¥ ;  this  defines  stationary  states  in  which  the 
phase  (and/or  time)  is  deliberately  left  unspecified.  Time  may  be 
pertinent  in  the  case  of  an  electron  (when  one  could  substitute  the 
Fourier  transformation  of  the  moving  charge  in  real  time),  but  this 
pertinence  does  not  extend  to  a  neutron  moving  in  a  potential  well. 
Consequently,  it  would  appear  that  the  theorists  have  not  fully 
examined  the  temporal  ramifications  of  postulating  that  particles 
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obey  such  wave  equations  as  the  Klein-Gordon  equation 


( □  2  +  m2)  V  ( x ,  t)  =  0  , 


(27) 


where  the  further  postulated  periodic  nature  of  the  temporal  para 
meter  implies  phase  not  time. 


THE  OPTICAL-MECHANICAL  ANALOGY 

In  the  previous  problems,  the  analogy  between  optical  and 
mechanical  solutions  is  most  striking.  Thus,  the  Hamiltonian  for 
the  mechanical  particle  in  spherical  coordinates, 

2^  (pr2  +  P#2  +  P(j)2  )  +  V  =  H  ,  on  substitution  of  the  operator, 

p  ->  -  i  —  (fi  =  1),  becomes,  in  =  E 'P 
o  q 


+ 


1 

r2  sin2  0 


-  2  m  E 


¥  =  0  , 


(28) 


closely  analogous  to  the  isotropic  ray  Hamiltonian 


o 


2 

r 


0  1  0  -  oq2  -  n2  (r,  0,< D)  =  0  . 

t2  sin^  $ 


(29) 


In  addition  to  delineating  the  role  of  spatial  isotropy  in 
quantum  theory  the  example  suggests  an  expanded  analogy  as,  for 
example,  in  Table  I. 
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TABLE  I 


Optical  Mechanical  Analogy 


Mechanical  Systems  (holonomic) 

Optical  Systems 

Canonical  conjugates 

Canonical  conjugates 

generalized  position  .  .  .  qk 
generalized  momentum  .  .  .  pk 

<9L 

k  ■ 

position  .  .  . 

“slowness”  •  •  •  cr^ 

“i- '  (i) 

Time  .  .  .  t 

Temporal  parameter  .  .  .  u 

phase  t  =  fa,  du 

J  K 

,i».  J  (n  *  „|i)  (-£-) 

Hamiltonian  .  .  .  H 

Ray  Hamiltonian  .  .  .  fi 

H  =  H  (pk  qk  ,  E  ,  t) 

fl  =  fi  (xk  ,  ak  ,  co  ,  t) 

for  time  stationary  systems 
dH 

ir  =  0 

for  time  stationary  media 
dfl 

*d r= 0 

Canonical  equations 

an  .  .  an 

a  pk  qk  ’  <9  qk  ~  Pk 

Canonical  equations 

an  .  an 

a<rk  xk-  axk  -  ~ 

Lagrangian  L  (qk  ,  qk)  =  T  -  V 

S^Ldt  =  0 

Indicatrix  v  (xk  ,  ak)  =  1 

sf(l)  ds  .  0 

Generating  function  of  the  Hamiltonian 

Generating  function  of  the  ray  Hamiltonian 
(Figuratrix) 

H  =  Pk  %  -  L 

0  V  CTi 

n  =  773 - 7 - r  -1  =  0 

*  & 3/  ® 9) 

from 

Atmospheric  Processes , 

Nawrocki  &  Papa 
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There  is  some  ambiguity  in  selecting  the  analogous  terms.  This 
may  come  about  as  a  result  of  over-defining  the  optical  system  by  the 
relation  E  =  -fuD.  Thus,  the  equations  of  motion  involving  the  two 
mechanical  canonical  sets  q^  and  H,  t^  can  actually  yield,  in  the 

optical  case,  only  a  unique  ray  trajectory.  Pertinent  to  this  am¬ 
biguity,  Landau  and  Lipshitz  apparently  prefer  to  use  as  the  canonical 
relations 


ki 


d  oj 


dco 


(30) 


where 
•— ¥ 

x  is  position 

— ♦  *-■> 

k  is  the  wave  vector  (k  n). 

— > 

Clearly,  the  wave  vector,  k  ,  and  the  frequency,  id  ,  transform  like 
momentum  =  ii k )  and  energy  (E  =  /ftuj) ,  respectively.  Using 

this  system  Landau  states 

"However,  it  (the  least  action  principle)  cannot  be 
written  as  6j  L  dt  =  0,  since  it  turns  out  to  be 
impossible  to  introduce,  for  rays,  a  function  analogous 
to  the  Lagrangian  of  a  particle  (i.e., 

L  =  p  •  —  H  ->  k  •  _  ty  =  0),  since  m  =  ck  in 

(9p  <?k 

vo lume . " 

While  this  is  true  in  the  sense  specified  by  Landau  ,  Synge 

[Ref.  6]  has  shown  that  the  optical  Lagrangian  can  be  formulated 

provided  that  a  slightly  different  procedure  is  followed.  He  starts 

dH 

with  H  defines  qp  =  —  —  »  solves  for  pp,  and  defines  L  as 
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L  (%  ■  %  •  *)  =  Pp  qP  -  H  (Pk  •  %  •  Q 
Thus,  Hamilton's  principle  for  the  photon  becomes 

8  j  (ffjj  dx^  +  Hdt  -  Hdt)  =  8  J  ak  dx^  =  0  , 

or  essentially,  6j  =  0,  i.e.,  Fermat's  Principle. 

The  complete  tabulation  of  the  analogy  as  specified  by  Synge 
where  c  =  h  =  1  ,  is  as  shown  in  Table  II. 


(31) 


26 


TABLE  II 

Optical  Mechanical  Analogy  as  Specified  by  Synge 


Term 

Optical 

Mechanical 

Coordinates 

xp 

qp 

Time 

x4  (=  it) 

t 

Momentum 

0P 

pp 

Energy 

aA  =  iH 

H 

Hamiltonian 

a  =  C74  —  iH  (dp  1  =  0 

H  =  H  (Pp  -  %  -  9 

Equations  of  motion 

idXp  d  il/d  Gp  ^ 

dx^  d^l/do^  dOp 

(?H 

V  "  dpp 

id 0^  d  S  x p  ^  j_j 

dx^  1  (?f 'l/doq  <9  Xp 

dH 

PP  '  <?qp 

Hamilton 

Variational 

Principle 

8  j° k  dxk=0 

8  J (P p  qp  -  H)  dt  =  0 

d  df  _  di 
ds  da k  <?xk 

(  Euler-Lagrange) 
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By  introducing  the  frequency  u>  as  the  fourth  member  in  the 
s lowness- frequency  tetrad,  Synge  has  demonstrated  that,  for  photon 
dynamics,  an  essentially  Loren tz -covariant  formalism  can  be  developed 
to  solve  such  problems  as 

"1.  A  set  of  unaccelerated  media,  separated  by  vacuum,  the 
optical  properties  of  each  medium  (when  at  rest)  being 
known,  the  media  may  be  heterogeneous  and  anisotropic. 

2.  As  above,  but  without  the  vacuum  separation,  the  media 
sliding  by  one  another  in  contact. 

3.  Fluids  in  accelerated  motion,  assuming  them  optically 
isotropic  in  the  local  or  rest  frame." 

While  it  is  acknowledged  that  optical  science  is  indebted  to  Synge's 

work  and  there  is  great  admiration  for  his  four-dimensional  construct, 

application  of  the  latter  is  limited  by  the  dynamic  character  of 

optical  media.  In  other  words,  all  ponderable  media  are  composed  of 

dynamic  mass  and  charge- bearing  corpuscles  whose  collective  behavior 

determines  the  optics.  As  a  result  of  these  mechanical  aspects 

strains,  slippage,  and  ruptures  occur  which  cannot  readily  be  included 

in  the  general  covariant  manner  appropriate  to  acceleration. 
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SECTION  III 


SCATTER  THEORY  OF  MICROWAVE  PROPAGATION 

STATISTICAL  OPTICS 

The  idea  of  developing  a  statistical  theory  of  rays  originated 
with  the  phenomenon  of  beyond- the-horizon  propagation  of  microwaves. 
Indeed,  most  of  this  section  was  written  in  1954  as  the  result  of  a 
suggestion  by  T.  F.  Rogers  to  seek  a  solution  (alternative  to  the 
Booker-Gordon  [Ref.  7]  solution)  that  would  provide  an  everpresent 
transmitted-signal  level  beyond  that  predicted  by  classical  propaga¬ 
tion  theory.  This  was  achieved,  but  the  theory  lacked  the  sophisti¬ 
cation  necessary  to  provide  for  such  important  parameters  as  bandwidth, 
fading  rates,  and  frequency  dependence. 

The  subject  of  this  paper  is  Hamilton  optics,  and  a  thorough 
review  of  the  scatter  theory  herein  would  not  be  justified.  Such  re¬ 
views  have  been  made  by  Staras  and  Wheelon  [Ref.  8],  Staras  [Ref.  9]  , 
Shkarofsky  [Ref.  10]  and  Vysokovskie  [Ref.  ll]  .  However,  existing 
scatter  theory  has  marked  capabilities  and  it  would  be  a  misrepresen¬ 
tation  to  dwell  only  on  its  failures;  consequently,  a  very  brief  re¬ 
view  of  the  salient  features  is  in  order. 

As  with  most  controversial  areas,  the  experimental  evidence 
varies  to  some  extent  with  the  investigator.  Again,  there  is  no 
attempt  to  weight  conflicting  numbers,  and  the  data  presented  must 
reflect  only  casual  acquaintance  with  this  area. 
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1.  Distance  Dependence 


The  field  intensities  as  a  function  of  distance  are  given  in 

Figure  2,  as  determined  by  T.  F. Rogers  et  al  [Ref.  12],  The  straight 

line  reveals  the  author's  opinion  that  the  mean  propagation  loss  is 

■*  les 

indeed  exponential  (e  ) ,  but  other  authors  are  just  as  convinced 
that  an  r  n  law  is  revealed  by  the  same  data. 

Ninety  percent  of  the  time,  the  median  field  strength  is 
approximately  57  db  below  free  space  at  100  miles,  with  a  further 
loss  of  12  db/100  miles.  For  one  percent  of  the  time,  the  attenuation 
relative  to  free  space  is  85  db  at  100  miles  and  again  increases  at 
0 . 1  db  per  mile. 


Figure  2.  Field  Intensities  for  Tropospheric  Scatter  as  a  Function 
of  Distance.  (after  T.  F.  Rogers  et  al) 
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2.  Frequency  Dependence 


Some  estimates  are  given  in  Table  III  but  the  author  adopts 
Bolgiano's  [Ref.  13]  estimate  that  the  frequency  dependence  varies  in 
time  (and  of  course  as  the  meteorological  character  of  the  air  mass) 
as  indicated  by  Figures  3  and  4. 

Bolgiano  also  indicates  (Figure  5)  a  relation  between  wave¬ 
length  dependence  and  dynamic  stability.  To  exclude  this  would  again 
be  deliberate  misrepresentation,  but,  at  this  juncture,  the  author's 
theory  is  not  sufficiently  sophisticated  to  make  quantitative  use 
of  the  dynamic  character  of  the  air  mass. 

3 .  Fading  Rate  and  Bandwidth 

As  indicated  in  Table  2,  there  are  two  characteristic  fading 
rates:  a  fast  component,  0.1  to  10  c/s,  which  shall  be  related  to 

the  time  required  to  alter  n  fluctuations  (with  scale  of  the  order 
of  millimeters)  by  collisional  processes;  and  a  slow  component,  10 
minutes  to  1  hour,  which  shall  be  related  to  dynamic  meteorologic 
processes  (windblown  clouds,  turbulence,  etc.).  The  bandwidth  will 
not  be  considered  at  this  initial  stage. 

Existing  Theory  of  TSM  (Tropospheric  Scatter  Mode) 

To  date,  theories  proposed  to  explain  the  TSM  phenomenon  have 
been  based  on  one  of  two  basic  concepts. 

1.  The  received  signal  is  due  to  scattering  from  turbulence- 
produced  inhomogeneities  in  that  medium  common  to  the  transmitting  and 


receiving  beams. 
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Figure  3.  Distribution  of  Wavelength  Dependence 

Differences  in  measured  hourly  median  basic  transmission  loss,  L, 
scaled  frequency  experiment,  417.05  Mc/s  versus  2290  Mc/s.  Round 
Hill-Crawf ord ' s  Hill  February  11,  1957  -  July  11,  1957,  241  cases 
(after  Chisholm,  Roche,  and  Jones,  1957). 
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Figure  4.  Seasonal  Variation  in  Distribution  of  Wavelength  Dependence 


Differences  in  measured  hourly  median  basic  transmission  loss,  AL, 
scaled  frequency  experiment,  417.05  Mc/s  versus  2290  Mc/s.  Round 
H ill--Crawf ord ' s  Hill:  Winter  --  February  +  March  1957,  76  cases 

Summer  --  June  +  July,  80  cases 
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Figure  5.  Variation  of  Wavelength  Dependence  with  Dynamic  Stability 
of  the  Atmosphere  in  the  Scattering  Volume 
r  =  0.80 


34 


Some  Characteristics  of  Scatter  Propagation 
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2.  The  received  signal  is  due  to  coherent  partial  reflection 
(normal  mode)  from  a  refractive  index  stratified  under  gravitational 
interaction. 

While  it  is  believed  that  both  of  these  concepts  do  pertain  at 
times,  in  varying  degree,  the  scatter  concept  will  be  given  more 
weight  here  since  it  dominates  current  thinking.  There  is  some 
justification  for  this  circumstance,  since  the  scatter  concept  has 
been  developed  to  a  far  more  sophisticated  stage  and  can  provide 
such  system  paramerters  as  fading  rate  and  bandwidth.  Of  course, the 
values  to  be  assigned  to  these  parameters  are  ad  hoc,  but,  in 
principle,  the  theory  does  have  this  extended  capability. 


Figure  6.  Basic  Scatter  Geometry 
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Legend: 


R^  =  distance  from  scattering  volume 
to  the  transmitter 

=  distance  from  scattering  volume 
to  the  receiver 

E.  and  E  are  wave  amplitudes  of 
me  s  r 

incident  and  scattered  waves 

V  includes  all  volume  elements  which 
contribute  to  the  scatter 

\  =  wavelength  of  radiation 

0  =  angle  through  which  the  scattering 
proceeds 

p  =  angle  between  the  direction  of 
propagation  of  scattered  wave  and  E^^ 

Figure  6  shows  the  scatter  geometry  that  determines  the  boundary 
conditons  for  solutions  of  the  wave  equation.  For  a  medium  of 
dielectric  constant 

*  (r,  t)  =  <0  +  Af  ( r ,  t)  ,  <32) 
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Balser  [Ref.  15]  has  shown  that  the  first-order  wave  equation  is 


V2E; 


+  k  2 
j  +  Ko 


(33) 


where 


w2  fo 


€ 


and  indicates  a  rectilinear  coordinate. 

The  solution  of  this  quasi-static  wave  equation,  subject  to  the 
appropriate  boundary  conditions,  gives  the  following  tabulated 
parameters . 


Scattered  Field  Strength 
^  2 


Es(R2)  4W^R,R 


oKl  k2 


J .... 


sin  p  Af  (r,  t)  exp  (ik  •  r)  d'3  r 


(34) 


where  k  =  -  k^ 

Correlation  Function 


C(£)  =  Af  frj)  A*  fr2)  /  (A*)2 


(35) 


where  a  bar  over  a  term  denotes  the  mean  value. 


V  (A  (y 


\  Ae(r)Af(r" 


+  R) d3  r 


2  773  ( A  f  2 f  2) 


f  -»  ik  • 
J  S(k)e 


R  o 

d3k 


Spectrum  Function 


s  (k)  =  J”  e* k  ‘  R  C  ( R)  Af2/ 102  d3R 
V 


(36) 


Spectrum  Density  Function 

d>  (k)  =  S  (k)  k2/2rr2 


(37) 


Scattering  Cross  Section 


a{d,  A,  p)  =  (R22/V) 


(R2) 


"incid 


=  t r2  A2  sin2  p 


J 


R  Af  -  j3 

K  ~2~  C  ( R)  d3  R 


(38) 


Scattered  Power 


PR/PT  =  77~2 


k/ 


gt  gr 


16  ^  y  Rj  R2 


2V2  ^l'A>d  ' 


(39) 
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where  K=  (l+  |  P|)^ 

P  =  the  reflection  coefficient  of  the  ground  at  grazing  incidence 
=  normalized  power  gain  of  transmitting  antenna 
=  normalized  power  gain  of  receiving  antenna 

While  acknowledging  the  great  flexibility  and  value  of  the 
Booker-Gordon  theory,  its  ad  hoc  nature  is  of  some  significance 
and  worth  stating  here.  When  the  theory  was  first  proposed,  the 
scale  of  the  typical  fluctuation  as  measured  at  the  University  of 
Texas  was  of  the  order  of  58  meters.  This  is  apparently  a  derived 
figure  somewhat  dependent  upon  the  correlation  function  assumed. 
Insertion  of  the  measured  Ae/e  for  inhomogeneities  of  this  size, 
together  with  an  anticipated  Gaussian  correlation  function,  led  to 
a  TSM  signal  which  was  orders  of  magnitude  smaller  than  that  ob¬ 
served.  Consequently,  Booker  and  Gordon  decided  to  select  an  ex¬ 
ponential  correlation  function  (this  yields  a  higher  TSM  signal) 
and  to  postulate  a  heirarchy  of  smaller  fluctuations  (micro turbulence) 
supported  by  the  larger  measured  scale  values.  Opponents  hastened 
to  point  out  the  heuristic  difficulties  with  the  exponential  corre¬ 
lation  function  (the  existence  of  a  cusp)  and  further  noted  the 
convenience  of  the  microturbulence.  The  latter  could  account  for 
the  measured  field  intensity,  and,  at  the  same  time,  was  of  such  a 
nature  that  its  existence  could  not  be  checked  experimentally. 
Nevertheless,  micro turbulence  was  a  step  in  the  right  direction. 
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When  it  was  also  found  that  obviously  quiet  atmospheres  gave 
significantly  higher  signals  (even  the  transmission  of  UHF- video  over 
path  lengths  of  200  miles),  a  further  ad  hoc  adjustment  was  required. 
This  was  provided  by  Bolgiano  (see  Figure  7)  who  varied  the  spectral 
shape  and  the  point  of  viscous  cutoff  T|  in  the  refractive  index 
fluctuation  spectrum. 

This  ad  hoc  character  of  the  theory  is  in  itself  no  great 
flaw.  However,  it  is  thought  that  Booker  and  Gordon  have  over¬ 
estimated  the  contribution  of  single  scatter  in  the  common  volume. 
Consequently,  their  concept  has  some  elements  quite  unreconcilable 
with  the  dynamics  of  the  atmosphere.  For  example,  it  led  quite 
naturally  to  Belgian1  s  suggestion  that  T|  varied  temporally,  lying 
in  the  meter  region  much  of  the  time.  This  is  contradicted  by  the 
meteorologists  who  estimate  T)  to  be  of  the  order  of  millimeters, 
a  result  motivated  by  our  own  considerations  of  tying  the  ever¬ 
present  microfluctuations  to  the  variations  in  the  Maxwell-Boltzmann 
distribution.  It  is  not  entirely  fortuitous  that  this  same  size 
microscale  is  motivated  by  the  plasma  physicists  as  that  (basically 
the  Debye  length)  appropriate  to  the  ionosphere.  As  a  possible 
motivation  for  this  it  is  suggested  that,  in  a  neutral  corpuscular 
medium,  molecules  comprising  a  microfluctuation  interact  through 
matter  waves  until  a  dimension  of  the  order  of  a  millimeter  is 
reached.  In  a  plasma  such  as  the  ionosphere,  interdependence 
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Figure  7.  Spectrum  of  Mean  Square  Refractive  Index  Fluctuations 


through  coulombic  interaction  (virtual  electromagnetic  waves  if  you 
like)  holds  until  a  similar  value  is  attained. 

The  nonphysical  nature  of  Booker's  view  on  turbulence  has  led  to 
the  suggestion  of  other  correlation  functions  (see  Table  IV).  However, 
each  of  these  models  preserves  the  concept  of  the  common  volume.  The 
next  section  shows  that  the  observed  everpresent  exponential  decay  is 
consistent  with  a  concept  requiring  neither  turbulence  nor  stratifi¬ 
cation,  but  requiring  the  active  volume  to  include  the  entire  region 
from  the  optical  horizon  to  the  receiver. 

Hamilton  Statistical  Optics 

An  attempt  to  devise  statistical  optics  was  influenced  by  the 
optical-mechanical  analogy  and  an  intuitive  appreciation  of  the  propa¬ 
gation  process.  However,  analysis  of  the  scatter  theory  of  propagation 
led  the  author  to  take  exception  to  the  concept  that  only  the  volume 
common  to  the  unscattered  transmitting  and  receiving  beam  is  impor¬ 
tant.  This  may  be  acceptable  for  a  localized  turbulent  medium  such 
as  a  rocket  plume  or,  to  a  lesser  degree,  ionspheric  scatter,  but  for 
the  case  of  tropospheric  scatter  both  the  transmitter  and  receiver  are 
enclosed  in  the  optical  instrument. 

Intuitively,  then,  it  is  the  author's  belief  that: 

(1)  The  received  signal  is  the  sum  total  of  rays 
that  extend  from  T  to  R. 

(2)  Each  ray  follows  a  path  as  described  by 

Hamilton's  equations. 
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TABLE  IV 


Correlation  Functions  and  Spectrum  of  Turbulence 
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(3)  There  are  no  discontinuities  in  the  medium. 

The  refractive  index  n  is  a  continuous 
function  everywhere.  The  atomistic  nature 
of  a  gas  provides  for  finite  incremental 
changes  in  index,  but  these  are  replaced 
by  a  fictitious  fluid  of  continuous  index 
which  provides  the  same  ray  trajectory. 

In  constructing  a  statistical  theory  of  rays,  the  form  of  the 
ray  Hamiltonian  (figuratrix  Q(xk,cjk,t)  =  0)  will  be  considered 
to  be  homogeneous  of  degree  1,  so  that  the  temporal  parameter  u 
is  the  phase. 

Introducing  phase  space  $  and  density  function  D,  the  total 
number  of  rays  is 

N  =  J  Dd<I> , 

where  d  =  D  (xk  ,  ak ,  t)  . 

If,  further,  D  =  D(Q),  then,  on  insertion  of  the  canonical 

equations 

aft 

d  xk  ° k 

a  ft 

a^  =  *k'  (40) 
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there  results 
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Consequently,  if  the  density  function  in  phase  space  is  a  function 
only  of  the  ray  Hamiltonian,  then  D  is  stationary,  and  Liouville's 
theorem  pertains. 

This  time  (or  rather  phase)  stationary  ensemble  of  rays  is 
assumed  in  the  motivation  of  a  propagation  law.  The  time- stationary 
mechanical  ensemble  is  the  Boltzmann  distribution 


P  (H)  =  C  e~H/KT  =  C  e” E/KT  , 


(42) 


2 

where  E  =  p  /2m  are  the  eigenfunctions  (here  continuous)  of  the 
Hamiltonian  H.  In  pursuing  the  analogy,  there  is  no  difficulty  in 
replacing  derivatives  of  H  by  the  corresponding  derivatives  of  tl, 
but  the  direct  substitution  of  U  for  H  produces  difficulty  since 
it  has  been  assumed  that  the  relation  ti(xk,(Tk,t)  =  0  always  holds. 
However,  the  concept  of  a  probability  function  P(Q)  =  P(xk,aR)  in¬ 
tuitively  exists,  and  it  is,  uniquely  determined  by  the  motion  of 
Par^^-C^es  comprising  the  Boltzmann  ensemble.  Since  only  the  time- 
stationary  case  is  being  considered,  no  turbulence  or  spectrum  of 
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eddies  are  present.  Only  the  time- stationary  particle  distribution 
is  involved. 

Assuming  a  very  large  number  of  molecules  N  enclosed  in  a 

fixed  volume  V,  let  v  be  a  fixed  subvolume  large  enough  to  hold 

a  large  number  ^m(o  «  m  «  N)J  of  molecules.  The  mean  value  of  m 

in  v  is  m  =  pN  and  the  actual  value  is  m  =  m  +  y.  Then, 
o  o 

by  the  Bernoulli  theorem, 


fN«  -  [2„n,0(l-m0/N,l  14  exp  {-  ^  (f _  %/H)  }  • 
If  N/mo  =  oo  ,  the  customary  relation 

foe  (y)  =  (2*mor'/2  exp  (  -  ) 


(43) 


(44) 


results}  but,  since  N  is  finite,  there  are  two  contributions  with 
standard  deviations 


From  these  particle  distributions,  the  phase  refractive  index 
f iguratrix  will  exhibit  mean  values  (termed  eigenvalues  of  the 
f iguratrix ) 
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Figure  8.  Figuratrices 

and  deviations  Aa^^  which  depend  upon 

1.  The  total  number  in  the  local  region,  v  (Gaussian 

° » 

f luctuation) . 

2.  The  number  of  weighted  constituents  such  as  H90  molecules 
(Gaussian  fluctuation)  in  the  local  region,  v  . 

3.  The  velocities  of  the  particles  (Boltzmann  distribution 

2 

“  mv  /  2KT 

of  energies  e  ,  again  Gaussian  distribution  of 

*  velocities) . 

The  sum  total  of  these  is  a  Gaussian  distribution  in  refractive  in¬ 
dex  eigenvalues  a^  . 

Gravitation  can  be  included  to  give  density  dependence  on 

altitude,  and  temperature  can  be  arbitrarily  imposed  according  to 
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the  customary  lapse  rates.  There  is  one  further  imposition  on  the 
fluctuations  in  phase:  the  phase  fluctuation  cannot  assume  such 
a  value  that  the  related  group  velocity  would  exceed  the  velocity 
of  light;  i.e., 


Vgroup 


c 

"group 
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n  +  o) 


d  n 
d  co 


<  c  . 
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Figure  9.  Illustration  of  Beyond -The-Horizon  Transmission 

Figure  9  shows  the  mean  path  for  an  air  mass  under  gravitation,  and 
tries  to  illustrate  how  the  Hamilton  rays  starting  at  T  actually 
fan  out,  each  ray  becoming  exponentially  less  probable  as  the  de¬ 
viation  from  the  mean  index  increases. 

If  it  is  assumed  that  the  rays  arriving  at  R  follow  a  path  of 
constant  curvature  (constant  deviation  Aa^)  ,  then  the  probability 
of  arriving  at  R  is  proportional  to 


49 


(47) 


i  =1 


This  is  exactly  the  expression  for  the  average  everpresent  signal. 


Idempotency  of  the  factors  apparently  implies  that  the  overall 

fluctuations  are  not  being  conserved.  Actually,  the  concept  of  a  con¬ 
stant  curvature  is  an  additional  postulate  because  it  is  clear  that 
Hamilton's  rays  exist  where  the  curvature  does  change.  Nevertheless, 
the  theoretical  view  predicts  that,  even  in  a  constant-temperature 
time- stationary  medium,  there  exists  a  great  many  small-scale  fluc¬ 
tuations.  Although  there  is  no  obvious  physical  mechanism  in  the 
neutral  gas  that  cap  play  the  role  of  the  Debye  length  in  ionized 
media,  the  mean  fluctuation  is  measured  in  millimeters  and  the  re¬ 
fractive  index  deviation  in  the  Gaussian  distribution  is  basically 
proportional  to  the  density,  (Figure  10).  Consequently,  the  idea  of  an 
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Figure  10.  Two  Representations  of  Exponential  Loss  in  Scattering 
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average  curvature  for  the  of f -the-median  trajectory  is  logical 
enough,  since  no  gravitational  separation  in  the  troposphere  is 
assumed  and  the  lower  part  of  the  troposphere  is  not  properly  weighted. 


Figure  11.  Meteorological  Conditions  Exhibiting 
A  Varying  Degree  of  Coherence 

In  Figures  11  (a  through  d),  cloud  forms  exhibit  an  increasing 
degree  of  coherence.  Both  cohered  refraction  and  partial  reflection, 
Carroll-Ring  [Ref.  16^,  and  turbulence,  Booker-Gordon,  can  then 
materially  alter  the  signal  characteristics.  In  spite  of  the  appar¬ 
ent  similarity  between  the  Hamilton  motivated  results  and  the 
scatter  theory  employing  a  Gaussian  correlation  function,  there  is 
a  basic  difference  due  to  the  assumed  size  of  the  scattering  volume. 
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When  an  isotropic  single  scatter  in  the  medium  common  to  the 
two  beams  is  assumed,  $  is  an  ignorable  coordinate  in  the  scattered 
amplitude.  Although  the  author  suggested  a  decade  ago  that  such  a 
^-dependence  should  be  expected,  the  interpretation  of  a  ^-dependence , 
if  found,  is  undoubtedly  ambiguous. 

There  are  of  course  more  meaningful  tests  of  the  validity  of 
the  intersecting- volume  model. 


1.  Crossed  beams  (Figure  12)  are  adjusted  to  intersect  a 
common  volume  V.  Existing  scatter  theories  would  indicate  a  high 
correlation  between  simultaneous  (in  the  Newtonian  sense)  signal 
intensities  and  fading  rates.  The  'Hamilton  ray  theory  would  predict 
the  low  degree  of  correlation,  approximately  equal  to  the  fraction, 

common  path _ . 

total  path  (~HR) 

2.  The  Hamilton  theory  states  that  the  rays  retain  phase  in¬ 
formation;  therefore,  this  phase  information  can  be  regained. 
Ideally,  for  the  common-volume  theories,  the  experiment  is  given 
schematically  in  Figure  13. 
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Figure  13.  Technique  for  Evidencing 

Phase  of  the  Optical  Instrument 

If  an  optical  frequency  (laser)  were  used,  the  resultant  of  the 
scattered  and  unscattered  signals  (arranged  to  be  of  the  same  mean 
amplitude)  could  be  photographed  and  detailed  information  could  be  de¬ 
rived  concerning  the  optical  instrument.  For  microwaves,  conceivably 
the  screen  could  be  replaced  by  a  large  matrix  of  silicon  diodes,  but 
unfortunately,  the  signal  levels  (70  to  100  db  below  free-space  levels) 
are  too  low  to  be  readily  detectable  by  such  a  means. 

The  Hamilton-ray  approach,  even  without  extending  its  sophisti¬ 
cation,  has  been  of  value  in  evaluating  the  various  scatter  theories. 
The  investigation  has  motivated  several  changes  in  the  Booker-Gordon 
theory.  These  changes  are  described  below. 

First,  turbulence  is  not  necessary  to  provide  a  scatter  signal. 
Even  a  constant- temperature  homogeneous  gas  has  the  everpresent  micro¬ 
structure  (millimeter  fluctuations)  capable  of  providing  beyond- the- 
horizon  transmission. 
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Second,  all  peripheral  considerations  support  the  choice  of  a 
small-scale  Gaussian  correlation  over  the  expotential  selected  by  Broker 
and  Gordon.  For  contributions  of  turbulence  to  the  scatter  signal, 
the  Obukhov-Silverman  j^Ref .  17 j  construct  is  preferred  to  that  of 
B  ooker-Gordon . 

Third,  the  model  of  the  intersecting  volume,  while  possessing 

great  simplicity,  cannot  provide  for  the  observed  exponential  propa- 
— k  s 

gation  law,  e  .  Only  a  large  number  of  successive  small-angle 
refractions  or  scatterings  in  the  entire  volume  can  do  this,  especial¬ 
ly  those  between  the  optical  horizon  and  the  receiver. .  This  concept 
is  supported  by  consideration  of  the  inverse  scattering  problem  be¬ 
cause  the  reflected  signal  provides  for  a  unique  dielectric.  I 


Hamilton  Backseat ter 


So  far,  the  propagation  medium  has  been  treated  as  being  contin- 
tinuous,  and,  this  has  produced  a  picture  of  the  forward-scattering 
process.  If  the  process  is  extrapolated  to  backscatter,  the  inter¬ 
pretation  of  the  ray  making  a  180°  turn  does  not  appear  to  be  in¬ 
tuitively  a  very  satisfying  concept.  Consequently,  it  would  be  de¬ 
sirable  to  readmit  the  existence  of  discontinuities  where  the  Hamilton 
jump  conditions  Act^  &x^  =  0  provide  for  the  mean  ray.  The  back- 
scattered  signal  can  then  be  obtained  by  the  construct  developed  by 
the  New  York  University  group.  The  radar  return  is  the  sum  of  extrema 
TT  obtained  by  wrapping  Hamilton  rays  about  each  discontinuity. 
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Figure  14.  Ray  Configuration  for  Backseat ter 


In  figure  14,  an  inverse- square  attenuation  is  taken  from  T  to  the 
tangential  points  A  and  B.  From  A  to  B  (and  B  to  A)  an  exponential 
loss  is  imposed.  This  mode  will  not  be  discussed  in  detail,  since 
plasma  physics  (interactions  of  waves  and  particles)  has  already  given 
a  satisfactory  solution  for  the  most  interesting  application,  the 
ionosphere.  However,  there  has  been  an  unwarranted  inclination  among 
radio  propagationists  to  give  credence  to  Gordon's  non-physical  con¬ 
cept  of  individual  electron  backs cattering  macroscopic  waves 
(X  >  XD  >  r  )  with  a  cross  section  consistent  with  a  unit  albedo 
spherical  reflector  of  classical  radius  rg  .  Such  a  concept  is 
fraught  with  difficulties  for  the  interdisciplinary  physicist. 

*  X  D  is  the  Debye  length. 
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Without  elaboration,  the  author  merely  reasserts  a  position  taken  in 
1958  that,  in  this  radar  case,  everpresent  Gaussian  fluctuations  in 
electron  density  provide  the  scattering  mechanism.  As  the  Hamilton 
trajectory  theory  indicates,  a  further  decrease  in  the  wavelength 
until  X  >  X  does  not  provide  for  a  radical  change  in  the  propa¬ 
gation  mechanism.  It  is  somewhat  paradoxical  that  Gordon  abandons 
the  Booker-Gordon  model  in  that  example  where  it  is  least  inappro¬ 
priate.  Indeed,  the  observed  radar  cross  section  can  provide  a  test 
of  the  appropriateness  of  their  model  of  intersecting  beams  in  the 
forward-scatter  case. 


i 
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SECTION  IV 


QUANTUM  OPTICS 


Quantum  theory  is  essentially  a  discipline  based  on  the  DeBroglie 
concept  that  all  microscopic  extracules  exhibit  both  corpuscular  and 
wave  characteristics.  Consequently,  the  quantum  theory  of  both  par¬ 
ticles  and  waves  has  a  common  element.  Indicative  of  this,  Dirac 


[Ref.  18^  has  shown  that  the  Schrodinger  representation  for  a  set  of 


equivalent  particles  is  identical  to  the  quantum  field  representation 
of  a  system  of  waves.  In  at  least  one  direction  these  ideas  have  been 
elegantly  developed.  In  the  inverse- scattering  problems  (determining 
the  character  of  a  dielectric  from  the  reflection  coefficient),  Kay 
and  Moses  [Ref.  141  employ  mapping  of  the  electromagnetic  wave  equation. 


(48) 


into  the  unidimensional  wave  equation 


+  V  (x)  ¥  (x)  =  P2  V  (X)  . 


(49) 


The  Gelf and- Levitan  algorithm  is  used  to  compute  the  potential  function 
and  the  wave  function  from  the  scattering  coefficient,  and  mapping 
provides  the  solution  of  the  electromagnetic  problem. 

In  theory,  the  backscatter  problem  discussed  in  the  previous 
section  could  be  treated  by  starting  with  the  observed  signal  and 
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developing  the  character  of  the  dielectric,  but  the  necessary  con¬ 
dition  of  a  continuous  dielectric  constant  may  be  overly  restrictive. 
Fortunately,  the  approach  that  motivates  a  dielectric  from  the  point 
of  view  of  the  plasma  physicist  and  then  derives  the  backscattered 
signal  has  provided  an  excellent  solution. 

Because  of  these  physical  considerations,  the  synthesis  of 
quantum  optics  by  analogy  could  not  be  expected  to  lead  to  anything 
new.  Nevertheless,  on  cursory  inspection  of  the  canonical  nature  of 

x  and  a  it  might  be  thought  that  there  is  a  limit  to  the  ability 

K.  K. 

of  an  observer  to  know  simultaneously  both  the  position  and  slowness 
of  a  ray.  From  the  quantum  mechanical  relations 

d 

pk  =  -  ™  TI 

d% 

P^k  “  ‘IfPk  =  “  5fk 


Pk  =  +r  [H,  Pk]  , 


(50) 


the  following  could  immediately  be  written 

.  -fi  d 
“  1  e  «9xk 


°i  xk 
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(51) 


Considerations  of  dimensional  identity  demand  inclusion  of  a  factor 
with  the  character  of  the  inverse  of  energy,  while  the  physics  of  the 
problem  imply  uniquely  that  this  energy  must  be  that  of  the  ray, 

E  =  'fKO(N).  Consequently,  Planck's  constant  vanished  from  the 
commutation  relation 


°l  xk  -  xe°k  =  -  TZ7 


Usually,  the  argument  is  reversed  so  that  Hamilton's  relations  in 
optics  and  mechanics  lead  to  quantum  theory,  provided  that  the 
Planck  relation  E  =  is  assumed. 

For  purposes  of  this  discussion,  the  basic  characteristics  are 
revealed  in  a  comparison  of  the  uncertainty  principle  for  mechanics 


Apk  A£  =  hS£k  , 


(52) 


and  the  analogous  expression  for  optics 


(53) 


This  demonstrates  that,  in  a  sense,  quantum  theory  was  implicit  in 
the  classical  notion  that  the  energy  of  a  ray  of  known  frequency  was 
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spread  over  the  wavelength  X.  As  may  have  been  anticipated,  quantum 
optics  of  the  ray  type  is  then  an  essential  redundancy. 

This  section  concludes  with  several  pertinent  peripheral  con¬ 
siderations.  First,  in  the  concept  of  quantum  optics,  the  author  has 
attempted  to  push  the  discipline  past  the  limitation  that  A  must  be 
large  compared  to  atomic  dimensions.  Actually,  as  X  decreases, 
atomic  collisional  processes  such  as  photoionization  become  important. 
Second,  only  geometric  or  ray  optics  have  been  considered.  The  mater¬ 
ial  does  not  bear  at  all  on  the  subject  of  quantum  field  theory  where 
operators  can  be  introduced  for  the  electromagnetic  field  variables, 

E,  H,  B,  D  via  canonical  potential  functions  ¥,  TT .  Third,  the 
theory  is  linear;  ray  intensities  are  so  low  that  no  change  in  the 
refractive  index  is  introduced.  This  automatically  eliminates  such 
phenomena  as  high- intensity  rf  or  laser  beams  where  ionization 
(breakdown)  may  be  produced  even  though  -fiy;  is  not  sufficiently 
large  to  accomplish  ionization  directly  (intermediate  activated 
states  play  a  prominent  role). 

Finally,  the  conclusion  that,  in  essence,  quantum  mechanics 
combines  the  concepts  of  classical  mechanics  and  optics  is  neither 
entirely  trivial  nor  barren.  For  small  systems  where  the  limit  to 
observation  is  relevant,  Dirac's  [Ref.  18]  concept  of  "superimpos- 
ability  of  translational  states"  was  considered.  To  explain  the 
experimental  results  from  interference  in  double  slits,  interference 
in  split  beams,  and  polarized  photons  at  oblique  incidence  in 
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tourmaline  crystals,  Dirac  suggested  that  translational  states  (for 
example,  photons  in  transit  between  observable  interactions)  carry 
what  are  classically  considered  as  mutually  exclusive  parameters. 

In  the  split-beam  experiment,  for  example,  roughly  half  of  the  photon 
traverses-each  beam  and  the  photon  then  interferes  with  itself  on  the 
subsequent  superposition  of  the  two  beams. 


It  is  true  that  if  physics  is  defined  as  that  discipline 
appropriate  to  the  observable  world  (as  indeed  the  author  does)  then 
Dirac1 s  invention  is  unimpeachable;  for  if  a  device  is  introduced  to 
identify  which  beam  the  photon  is  really  (observably)  in,  then  the 
interference  effect  automatically  and  quite  logically  vanishes. 
However,  Dirac's  ingenious  construct  is  an  invention  hardly  to  be 
labeled  as  a  subjective  truth.  Personally,  the  author  finds  it 
disconcerting  to  admit  that,  even  beyond  the  ken  of  observation,  a 
photon  must  be  one  (a  corpuscle)  yet  cannot  be  one  (follows  two 
different  paths).  The  gymnastics  required  for  pursuit  of  quantum 
optics  now  bear  some  fruit  since  it  suggests  a  model  different  from 
that  of  Dirac's  superimposable  translational  states.  If  quantum 
mechanics  is  indeed  the  combined  optical-mechanical  discipline,  then 
all  wavicles  (e.g.,  photons  and  electrons)  are  associated  with 
Schrodinger  waves.  For  the  electron,  Heisenberg  [Ref.  25]  gives  as 
the  velocities  of  the  phase  and  group  waves  the  values 


E  p2/2m 

V/  =  —  =  — - 

0  p  p 


v/2)  , 
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where  E  is  the  kinetic  energy  of  the  electron. 

Thus,  the  phase  wave  associated  with  the  wavicle  has  the  velocity 

equal  to  one-half  the  group  wave  value  (v  is  of  course  the 

S 

velocity  of  the  corpuscle)  in  contradiction  to  the  author's  under 


kinetic  energy,  the  author  understands  function  of  the  total  energy, 


Since,  the  phase  wave  of  a 


>  c  . 


P 


wavicle  possessing  rest  mass  is  always  greater  than  the  characteristic 
velocity,  the  phase  wave  can  be  used  as  the  carrier  for  the  proba¬ 
bility  indicator  required  for  translational  states.  Each  time  the 
wavicle  interacts  with  matter,  the  information  is  impressed  on  the 
phase  wave.  The  wavicle  only  goes  one  path  (with  the  group  velocity 
v  )  but  in  the  subsequent  or  second  interaction,  the  new-distorted 
phase  wave  interacts  with  its  new  boundary  condition  to  produce  the 
phenomenon  of  self- interference.  If  a  third  boundary  condition  is 
placed  somewhere  in  the  path  between  incidents  1  and  2,  then  the 
phase  wave  is  further  distorted,  so  that  the  simple  phase  relation  of 
incidents  1  and  2  necessary  for  1,  2  self -interference  no  longer 
exists.  At  all  times,  a  photon  or  electron  (ignoring  coulombic 
forces)  interacts  with  itself  via  its  own  phase  wave.  The  result  is 
the  same  as  that  achieved  by  Dirac's  construct,  but  the  basic  self- 
contradiction  in  the  world-beyond-observation  is  removed. 
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Peripherally,  the  concept  of  a  phase  velocity  in  excess  of  the  char¬ 
acteristic  velocity  obviously  does  not  violate  relativity,  since  the 
observable  signal  is  carried  at  the  group  velocity  (still  the  velocity 
of  the  wavicle). 

These  considerations  support  the  notion  that  the  quantum- 
mechanical  t  can  be  either  phase  or  time,  the  proper  choice  being 
dependent  upon  the  specific  problem.  No  light  has  been  cast  upon 
the  question  of  the  appropriateness  of  spatial  isotropy  for  nucleons 
inside  the  nucleus. 
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